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ABSTRACT 

We study the reionization history of the Universe in cosmological models with non- 
Gaussian density fluctuations, taking them to have a renormalized \ 2 probability dis- 
tribution function parametrized by the number of degrees of freedom, v. We compute 
the ionization history using a simple semi-analytical model, considering various possi- 
bilities for the astrophysics of reionization. In all our models we require that reioniza- 
tion is completed prior to z = 6, as required by the measurement of the Gunn-Peterson 
optical depth from the spectra of high-redshift quasars. We confirm previous results 
demonstrating that such a non-Gaussian distribution leads to a slower reionization as 
compared to the Gaussian case. We further show that the recent WMAP three-year 
measurement of the optical depth due to electron scattering, r = 0.09 ± 0.03, weakly 
constrains the allowed deviations from Gaussianity on the small scales relevant to 
reionization if a constant spectral index is assumed. We also confirm the need for a 
significant suppression of star formation in mini-halos, which increases dramatically 
as we decrease v. 
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1 INTRODUCTION 

The reionization of the Universe is a direct result of the 
formation of luminous sources which produce the photons 
responsible (see Barkana & Loeb 2001 for a review). Given 
that the first stars appear in low-mass halos formed at high 
redshift, the reionization history of the Universe is expected 
to be a powerful probe of the amplitude and nature of den- 
sity fluctuations on small scales. 

The detection of Gunn-Peterson troughs (Gunn & Pe- 
terson 1965) in the absorption spectra of distant quasars 
suggests a late reionization for the Universe, at a redshift 
z ~ 6 (Becker et al. 2001; Fan et al. 2003; White et al. 2003; 
Fan et al. 2005; Gnedin & Fan 2006). The indications of a 
high optical depth from the first-year WMAP results (Kogut 
et al. 2003; Spergel et al. 2003), implying early reionization, 
led to a flurry of papers on quite complex ionization scenar- 
ios (e.g. Haiman & Holder 2003; Hui & Haiman 2003; Ciardi, 
Ferrara & White 2003; Cen 2003; Chiu, Fan & Ostriker 2003; 
Melchiorri et al. 2005). Simplicity has now largely been re- 
stored by the recent more precise third-year WMAP results 
(Page et al. 2006; Spergel et al. 2006), which indicate a much 
smaller optical depth, r = 0.09±0.03. This is consistent both 
with the quasar results and with a much simpler reioniza- 
tion history with a fast transition from a neutral to a ionized 
universe (see for example Choudhury & Ferrara 2006). Re- 



cent results by Haiman & Bryan (2006) suggest that this 
can only be achieved if star formation in high-redshift mini- 
halos has been significantly suppressed (see also Wyithe & 
Loeb 2006). 

In the present paper we perform a complementary 
study, investigating the dependence of the reionization his- 
tory of the Universe on the nature of the cosmological per- 
turbations, in the context of cosmologies permitted by the 
WMAP three-year results. We are motivated by the sugges- 
tions that a non-Gaussian contribution, for example from 
cosmic defects, might have been able to reconcile the high 
optical depth suggested by the WMAP first-year results 
and a low redshift of (complete) reionization implied by 
the quasar data (Avelino & Liddle 2004; Chen et al. 2003). 
With the new simpler ionization models, non-Gaussian mod- 
els should be better constrained. 

We compute the reionization history of the Universe 
using a simple semi-analytical model similar to that used 
by Haiman & Bryan (2006) (based on Haiman & Holder 
2003). We consider various possibilities for the astrophysics 
of reionization, subject to the constraint that reionization is 
completed prior to z — 6 as required by the measurement of 
the Gunn-Peterson optical depth from the spectra of bright 
quasars at high redshift. 
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2 THE MASS FRACTION 



3 THE POWER SPECTRUM 



We use the Press-Schechter approximation (Press & 
Schechter 1974) to compute the mass fraction, f(M), as- 
sociated with collapsed objects with mass larger than a 
given mass threshold M. This was originally proposed in the 
context of initial Gaussian density perturbations, and later 
generalized to accommodate non-Gaussian initial conditions 
(Lucchin & Matarrese 1988; Chiu, Ostriker & Strauss 1998). 
In both cases the mass fraction is assumed to be proportional 
to the fraction of space in which the linear density contrast, 
smoothed on the scale M, exceeds a given threshold <5 C : 



f(M) = Aj 



T(8)dS. 



(1) 



Here V(S) is the one-point probability distribution function 
(PDF) of the linear density contrast S, and Af is a constant 
which is computed by requiring that /(0) = 1 (Af = 2 in 
the case of Gaussian initial conditions). This generalization 
of the original Press-Schechter approximation has success- 
fully reproduced the results obtained from iV-body simu- 
lations with non-Gaussian initial conditions (Robinson & 
Baker 2000). 

A top-hat filter will be used to perform the smoothing 



W(kR) = 3 



sin(kR) cos(kR) 
~JkRf (kR) 2 



(2) 



where M — 47rp m 7? 3 /3 (here p m is the background mat- 
ter density). We take 5 C = 1.7, motivated by the spherical 
collapse model and ./V-body simulations. For spherical col- 
lapse 5 C is almost independent of the background cosmology 
(e.g. Eke, Cole & Frenk 1996). 

The evolution of the dispersion of the density field is 
given by 

g(f2 m ,fU) (7r(0) 



(3) 



where the suppression factor g(Slm, £^a) accounts for the de- 
pendence of the growth of density perturbations on the cos- 
mological parameters f2 m and J1a (Carroll, Press & Turner 
1992; Avelino & de Carvalho 1999). 

We shall assume that the non- Gaussian density contrast 
has a chi-squared probability distribution function (PDF) 
with v degrees of freedom, the PDF having been shifted so 
that its mean is zero (such a PDF becomes Gaussian when 
v — > oo). Hence, 



f(M) 



Q(is/2,S c /a R x + 
Q(v/2,v/2) 



(4) 



where Q(a, x) with a > is the incomplete gamma function 
defined by 



T(a,x) 



Q(a,x) = 
Here 
T(a,x) EE 

and T(a) = T(a, 0). When noowe have 
S c 



f(M) -> erfc 



V2 



(5) 



(6) 



(7) 



Historically, one of the main uncertainties in determining 
the reionization history of the Universe was the lack of an 
accurate determination of various cosmological parameters. 
However, these uncertainties have largely been removed by 
a growing body of precise cosmological data. Throughout, 
we shall adopt a cosmological model motivated by the three- 
year WMAP results (Spergel et al. 2006). As they find that 
WMAP alone does not require a running of the spectral 
index, we will take as our base cosmology their preferred 
model with a power-law initial spectrum. The parameters 
are a matter density = 0.24, dark energy density Q A = 
0.76, baryon density fig = 0.042, Hubble parameter h — 
0.73, normalization erg = 0.74, and perturbation spectral 
index n = 0.95. 

To determine the amplitude of perturbations on the 
short scales relevant to reionization, we use the transfer func- 
tion from Bardeen et al. (1986), so that the power-spectrum 
is given by 



P(jfc) oc k n 



ln(l + e q) 



-1/4 



where q = k/hV, [k] = Mpc~\ 
e= [2.34,3.89,16.1,5.46,6.71], 
and 



(8) 



(9) 



(10) 



is the shape parameter. The dispersion of the density field 
is given by 



2 

or 



-f 

Jo 



k 2 V{k)W 2 (kR)dk. 



(11) 



For Gaussian fluctuations the Press-Schechter formal- 
ism leads to an underestimate of the halo mass function at 
the high- mass end (Jenkins et al. 2001), while the presence of 
baryons leads to a suppression of power on small cosmolog- 
ical scales. These two opposite effects are of the same order 
of magnitude and consequently we shall ignore the baryon 
correction to the shape parameter (Sugiyama 1995). Also, 
it has been shown that small deviations from the predicted 
mass function could arise in particular in the limit of rare 
events (Avelino & Viana 2000; Inoue & Nagashima 2002). 
However, since we will fix the redshift at which reionization 
is to be completed by adjusting efficiency parameters, these 
uncertainties will not affect our results. 



4 THE REIONIZATION MODEL 

Following Haiman & Bryan (2006) (see also Haiman & 
Holder 2003) we will classify dark matter halos into three 
different categories according to their virial temperatures 

4 x 10 2 K < T < 10 4 K (Type II) , 
10 4 K < T^ 9 x 10 4 K (Type la) , 
~ T > 9 x 10 4 K (Type lb) . 

The total fraction of the mass in the Universe that is con- 
densed into Type II, la, and lb halos is given by 
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Fcoii,ii(z) = /(M„) - /(Mi a ), 
Fcoii,i a (2) = /(M Ia ) - /(Mi b ), 
and 

Fcoil,ib(z) = /(Mi b ), 



(12) 
(13) 

(14) 



where the halo masses Mn = Mii(Th, z), Mia, = A/i a (Ti a , z), 
and Mib = Mib (Tib, z) for the virial temperatures Tu = 



4 x 10 2 K, Ti 
from 



10 4 K and T m = 9 x 10 4 K are obtained 



M 
10 6 M f 



/ T \ 3/2 / l +£ \ 
U800AV V 21 y 



-3/2 



(15) 



These halos differ in their virial temperatures, and conse- 
quently in the cooling processes which are effective in each 
case. Ionizing sources can form in type II halos only in the 
neutral regions of the intergalactic medium and provided 
there is a sufficient abundance of H2 molecules. On the other 
hand, type I halos can cool and form ionizing sources inde- 
pendently of the H2 abundance. However, while type la halos 
can form new ionizing sources only in neutral regions, the 
formation of new ionizing sources in type lb halos can also 
occur in the ionized regions of the intergalactic medium. 

The recombination rate per unit of time and volume can 
be written as aBCHii(i-Hii) 2 , where qb is the recombination 
coefficient of neutral hydrogen to its excited states (= 2.6 x 
1CP 13 cm 3 s" 1 at T = 10 4 K), (nmi) is the mean number 
density of ionized hydrogen, and Chii = (^hii)/(^hii) 2 is 
the mean clumping factor of ionized gas. Assuming that Chii 
does not vary with redshift, the probability P(ti,t) that a 
ionizing photon produced at a time ti is still contributing to 
reionization at a time t > ti is given by 



P(ti,t) = exp 



(16) 



where i r = obCuuX pB(ti)tf /m p where pb is the baryon 
density, X is the hydrogen mass fraction and m p is proton 
mass. Alternatively we may adopt the simple power-law re- 
lation of the evolution of the dumpiness Chii advocated by 
Haiman & Bryan (2006) 



Chii(z) 



1 = 9 



1 + z 
7 



-0 



(17) 



which in the /3 = case reduces to Chii = 10. 

The HII filling factor Fuu(z) is then obtained from the 
collapsed gas fractions as 



/■1 nr.) - / dz' [ei^^iz') + [1 - Fun(z')] 



dz 



( Z ')+en^f^(z')]}p(z',z) 



(18) 



where ei and en are respectively the efficiencies (number of 
ionizing photons produced per proton in collapsed regions) 
of type I halos and type II halos (it being assumed that the 
efficiencies in type la and lb are equal). The second term 
on the right-hand side of equation 118H includes a factor of 
(1 — Thii), which explicitly takes into account the fact that 
new ionizing sources should appear in type II and type la 
halos only in regions which have not yet been ionized. The 
injection efficiency, e, of ionizing photons into the intergalac- 
tic medium can be parametrized as e = A 7 /»/ csc . Here /* is 
the fraction of baryons that turn into stars, iV 7 is the average 



number of ionizing photons produced per stellar proton, and 
/ csc is the fraction of these ionizing photons which escape 
into the intergalactic medium. 

Type II halos are expected to host very massive stars, 
which are efficient at generating ionizing photons. On the 
other hand, it is expected that only a small fraction of the 
available gas will be incorporated into stars. Type I halos 
are less efficient at producing ionizing photons, but that is 
compensated by a larger /». In this paper we will assume 
that the efficiencies do not depend on z and consider the 
values advocated by Haiman & Bryan (2006) for type II 
halos, that is A 7 = 80000, f, = 0.0025 and / csc = 1, so 
that en = 200 if no suppression of star formation in type 
two halos occurs (en = otherwise). We fix the efficiency 
of type I halos such that the Universe becomes fully ionized 
at z = 6.5. This means that e\ depends on v. However, 
we find that the efficiencies are lower only by a factor of 
approximately 2 in the v — 1 case compared to the v — > 00 
one. In the Gaussian case this approximately corresponds to 
the parameters chosen by Haiman & Bryan (2006), that is 
A 7 ~ 4000, /„ ~ 0.15 and / csc ~ 0.2 for type la and type lb 
halos so that ei ~ 120. 

The optical depth is calculated using 



t(z) = cctt / dt JXe(t') 



I.OSccttX f , dt 

dz — psFvax 
dz 



dz' (i±£ 



(nS.(i + z'r + n° A ) 1/2 



Fun. 



where 

r = — ~ 0.056 il B h. 

STvGm„ 



(19) 
(20) 

(21) 



The factor of 1.08 approximately accounts for the contribu- 
tion of Helium reionization. Here we assume that HII and 
Hell fractions are identical and that the Hell to Helll tran- 
sition takes place at 2 < 6. 



5 RESULTS 

In Figure 1 we plot the mass fraction in type II, type la and 
type lb halos, as a function of redshift z, for chi-squared 
distributions with different numbers of degrees of freedom 
(y — 1, 10, 100, 00). We clearly see that the evolution of mass 
fraction with redshift becomes becomes more gradual for 
smaller values of u, thus making reionization a slower pro- 
cess in the non-Gaussian case. This can be confirmed in 
Figure 2 where we plotted the evolution of the correspond- 
ing ionized fraction of hydrogen as a function of redshift z, 
for two possible ionization histories (see Figure 2). In each 
case we required that reionization is completed by z = 6.5. 
The solid line represents a reionization model with no con- 
tribution from type II halos (en = 0), while the dotted line 
represents a model where type II halos are the dominant 
ionization source at high redshift (en = 200). As expected, 
the slower evolution of the mass fraction in type la, type lb 
halos and type II halos with redshift at small u is reflected 
in the evolution of ionized fraction of hydrogen, which also 
becomes a slower function of redshift. 
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Figure 1. The evolution of the mass fraction, -Fcolli in type II, 
type la and type lb halos (solid, dashed and dotted lines respec- 
tively), as a function of redshift z, for chi-squared distributions 
with different numbers of degrees of freedom (y = 1, 10, 100, oo). 



Figure 3. The optical depth due to Thomson scattering as a 
function of redshift, z, for the models considered in Fig, l2l 
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Figure 2. The evolution of the ionization fraction of hydrogen, 
Fhii, as a function of redshift, z, for the mass fractions considered 
considered in Fig. HI and for two possible ionization histories, ej is 
chosen so that reionization is completed by z = 6.5. The solid and 
dotted lines represent reionization models with ejj = (negligible 
contribution from type II halos) and en = 200 (type II halos are 
the main ionization source at high redshift) respectively. Here we 
assumed that the clumping factor, Chii, is given by equation 1171 
with (3 = 2. 
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Figure 4. The evolution of the ionization fraction of hydrogen, 
Fhii, and the integrated optical depth due to Thomson scatter- 
ing, t, as a function of redshift, z, for a model with v = 1 (top) 
and a model with v = oo (bottom). Here we consider two possi- 
ble ionization histories with ejj = (negligible contribution from 
type II halos) and en = 200 (type II halos are the main ioniza- 
tion source at high redshift) and two different evolutions for the 
clumping factor parametrized by fi = 2 (solid line) and j3 = 
(dashed line). 



In Figure 3 we plot the optical depth due to Thomson 
scattering as a function of redshift, z. Here, we confirm the 
results of Haiman & Bryan (2006) showing that in order not 
to overproduce the optical depth, star formation in mini- 
halos appears to have been suppressed even if we assume 
Gaussian fluctuations. As we decrease v the need for this 
suppression increases dramatically. 

In Figure 4 we show the impact of two different 
parametrizations for the evolution of the clumping factor 
(/3 = 2 and /? = 0) on the reionization history of the Uni- 
verse. We see that the poor knowledge of the redshift depen- 
dence of the clumping factor, Chii(z), introduces significant 



uncertainties in the predicted reionization history, which be- 
come increasingly large at small v. 

Finally, in Figure 5 we plot the integrated optical depth 
(up to z = 60) due to Thomson scattering, r, as a function 
of the number of x 2 degrees of freedom considering different 
evolutions for the clumping factor parametrized by j3 = 2 
(solid line), [3 = 1 (dashed line) and = (dotted line) 
and a reionization model with a negligible contribution from 
type II halos (en = 0). The horizontal stripe represents the 
range of r allowed by the WMAP three-year data (at 68% 
confidence level) . We see that given the model uncertainties 
reionization alone cannot yet rule out any value of v with a 
significant degree of confidence. 
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Figure 5. The integrated optical depth (up to z = 60) due to 
Thomson scattering, r, as a function of the number of x 2 degrees 
of freedom considering different evolutions for the clumping factor 
parametrized by fi = 2 (solid line), fi = 1 (dashed line) and 
fi = (dotted line) and a reionization model with en = 0. The 
horizontal stripe represents the range of r allowed by the WMAP 
three-year data (at 68% confidence level). 



6 CONCLUSIONS 

We have investigated the possible impact of non-Gaussian 
density perturbations on the reionization history of the 
Universe in the context of cosmological models with non- 
Gaussian density fluctuations with a renormalized \ 2 prob- 
ability distribution function parametrized by the number, 
v, of degrees of freedom. We assumed a constant spectral 
index so that the amplitude of the density perturbations on 
the small scales relevant to reionization could be calibrated 
using the three-year WMAP results. We have confirmed the 
results of Haiman & Bryan (2006) suggesting a significant 
suppression of star formation in mini-halos. 

We have also shown that the reionization history of the 
Universe is able to constrain the non-Gaussian nature of the 
density perturbations. However, we have seen that reioniza- 
tion alone is not yet able to rule out any value of v (although 
very small values of v appear to be disfavoured) . Given that 
large deviations from a Gaussian probability distribution 
function on large cosmological scales are already excluded 
by the WMAP data, we conclude that unless the probability 
distribution function has a strong scale dependence the cos- 
mic reionization constraints are not competitive with those 
coming from the cosmic microwave background anisotropies. 
Still, it is important to bear in mind that the scales probed 
by reionization are much smaller than those probed by the 
WMAP (or even Planck), and consequently the consistency 
is reassuring, in particular taking into account recent claims 
(Mathis, Diego & Silk 2004) of hints for non-Gaussianity on 
scales much larger than those relevant for reionization. 

We should also emphasize that our assumption of a con- 
stant spectral index excludes an important class of models 
where the density field is the sum of Gaussian and non- 
Gaussian contributions with different power spectra. On 
large cosmological scales the non-Gaussian component is 
severely constrained and only very small deviations from 
Gaussianity are allowed. However, if the power spectrum of 



the Gaussian contribution is steeper than that of the non- 
Gaussian one, as happens in hybrid models with both in- 
flationary and defect perturbations, the non-Gaussian part 
may be the dominant component on small scales while being 
completely negligible on large cosmological scales. 

Avelino & Liddle (2004) have recently shown that the 
reionization history of the Universe leads to stringent con- 
straints on the energy scale of defects (see also Olum & 
Vilenkin 2006). However, since in that case the redshift at 
which the Universe becomes fully ionized no longer fixes the 
efficiencies of type I halos, the analysis is more dependent 
on the astrophysics of reionization. We shall return to this 
issue in a forthcoming publication. 
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